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Structure formation with cold dark matter (CDM) predicts halos with a central density cusp, 
which are observationally disfavored. If CDM particles have an annihilation cross section av ~ 
10" 29 (m/ GeV) cm , then annihilations will soften the cusps. We discuss plausible scenarios for 
avoiding the early Universe annihilation catastrophe that could result from such a large cross section. 
The predicted scaling of core density with halo mass depends upon the velocity dependence of av, and 
s-wave annihilation leads to a core density nearly independent of halo mass, which seems consistent 
with observations. 
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Introduction. The idea that the large-scale structure de- 
veloped by gravitational instability from initially small- 
amplitude, adiabatic and nearly scale-invariant fluctua- 
tions is compatible with a number of observables across 
a wide range of length scales (e.g., from the cosmic mi- 
crowave background anisotropy to the Lyman-a forest). 
Essential to this compatibility is the existence of cold 
dark matter: matter which is non-baryonic, has only very 
weak interactions with photons and baryons, and (prior 
to gravitational collapse) is cold. 

The greatest challenge to this otherwise success- 
ful scenario comes from the apparent discrepancy be- 
tween predicted dark-matter halo-density profiles and 
those inferred from observations. Simulations with non- 
interacting cold dark matter lead to halo density pro- 
files that are singular at the center |lj whereas ob- 
servations indicate uniform density cores. In this Letter 
we explore the possibility that the dark matter today 
has a large cross-section for annihilation which results in 
preferential destruction in high-density regions, softening 
halo cores ||. 

Detecting and determining the properties of the dark 
matter is a major goal of observational cosmology. If 
annihilations are indeed altering the properties of dark 
matter halos, then we have a new means of studying the 
dark matter. The interactions of the CDM particles de- 
termine both the magnitude and velocity dependence of 
the annihilation cross section. For example, for s-wave 
annihilation, o~a\v\ is independent of velocity and for p- 
wave annihilation, a a M is proportional to v 2 . These two 
different dependences result in different scaling relations 
between core density and halo velocity dispersion, which 
can be tested by current observations. 

As we show below, current data for high velocity dis- 
persion systems such as clusters of galaxies to low velocity 
dispersion systems such as galactic satellites are consis- 



tent with the same core density of about 1 GeV/ cm 3 (= 
0.026 M Q /pc 3 ) This scale invariance can be ex- 

plained by s-wave annihilation with a cross-section av ~ 
10 _29 (ro/ GeV) cm 2 , although future improvements in 
both the data and predictions will be necessary before 
such a statement can be made with confidence. As we 
shall discuss, the cosmological and astrophysical con- 
straints on annihilating CDM point to a candidate be- 
yond those currently favored (e.g., axion, neutralino). 
Halos of Annihilating Dark Matter. Numerical simula- 
tions of structure formation in the CDM scenario show 
that the dark matter halos which form with a wide range 
of masses are all well-fit with the so-called NFW Q form 
for the density profile. This form has p oc r -3 at large r 
and a cuspy inner region with p oc r~ a with a = 1. More 
recent higher resolution simulations || predict cusps that 
are even stronger, with a ~ 1.5. Nevertheless, in most of 
what follows, we use the NFW theory for simplicity. 
To be precise, the NFW profile is 

p(r = xr s ) = p s x _1 (l + x)~ 2 , (1) 

where the value of p s is determined by the mean den- 
sity of the Universe at the time the halo collapsed. In 
CDM theory, small objects collapse first, followed later 
by larger ones. Thus, there is an inverse relationship be- 
tween p s and halo size. In Fig. 1 we show this scaling 
relation with halo size represented by velocity dispersion 
for the halo, estimated as cr V i r = \J GM v \ x j2r w \ T . (The 
virial radius, r v ; r , is defined such that the mean den- 
sity inside the r v - n sphere is 200 times the present mean 
density of the Universe, and M v ; r is the mass contained 
within r v i r 0.) 

Annihilations will alter the halo profiles near the core 
where the density of the dark matter particles is the high- 
est. The annihilation rate (per particle) Y = n(cr|v|) de- 
pends on the velocity dispersion. We parameterize the 
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velocity dependence as T = (p/m)<7AV n (n = for s- 
wave; n = 2 for p-wave), where v is the velocity disper- 
sion and m is the CDM particle mass. 
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FIG. 1. The halo density at r = 0.1r s (where the cusp 
problem becomes prominent) in structures of different size 
according to NFW jl[ (solid curve). The objects are charac- 
terized by their virial velocity dispersion as indicated. An- 
nihilation lines (dashed curves), normalized to LSB galaxies, 
are shown for the cases n = and 2. Above the line, annihi- 
lations are very important (at r = 0.1r s ) and below the line 
they are unimportant. For n = 0, LSB and smaller objects 
have their cores softened significantly, while clusters do not, 
consistent with observations. For n — 2, clusters would be 
adversely affected. 

Fig. [I] also shows in a qualitative way how annihila- 
tions affect the core structure of different objects. The 
annihilation lines drawn show whether or not annihila- 
tions are important at a NFW halo radius of 0.1r s in 
different kinds of objects. (Note, since halo densities 
diverge, for any object annihilations become significant 
deep enough into the core). The annihilation lines are 
normalized to soften the cores of low surface-brightness 
(LSB) spiral galaxies. Because of how annihilations scale 
with velocity, for n = clusters remain unaffected at 
r > 0.1r s , while the cores of LSBs and smaller objects are 
dramatically softened. For n = 2, the opposite is true, 
which contradicts observations that indicate the NFW 
profile works well for clusters. We expect any n > 1 to 
be inconsistent with observations. The case of n = 0.5 
is interesting since the annihilation line runs parallel to 
the structure line (for cr v i r < 100 km/s), implying that 
all the systems will be smoothed off at the same value of 
r/r s . 

Model Building Constraints. For the annihilations to be 
effective in galaxy cores today, the annihilation rate must 
satisfy the (approximate) constraint: 

T ~ (p/plsb) (v/v LSB ) n H , (2) 

where the subscript LSB denotes the appropriate values 
for a typical LSB and Hq — lOO/i km s _1 is the present 



expansion rate of the Universe. Outside collapsed objects 
today, the density of CDM is much lower and annihila- 
tions will be unimportant for n > 0. The early Universe 
is another matter as densities were much higher, p cx T 3 , 
where T is the cosmic background radiation temperature. 

The figure of merit for the effectiveness of annihila- 
tions in the early Universe is measured by annihilation 
rate divided by the expansion rate: when T/H > I an- 
nihilations are effective (and vice versa). Assuming that 
the velocity dispersion of the CDM particles can be char- 
acterized by the background radiation temperature and 
normalizing the cross section to the desired value today, 
the temperature dependence of T/H is 

where T cq ~ I eV is the temperature at matter - radia- 
tion equality. There are three important things to note: 
(1) the large coefficient in front of this expression - an- 
nihilations in the early Universe are a significant consid- 
eration; (2) for n = — I , the effectiveness of annihilations 
is epoch independent and disastrous; and (3) for n > — 1 
annihilations were more important in the past. 

Observational data suggest that if halos are made of 
annihilating CDM particles, their annihilation cross sec- 
tion is characterized by n, < I. Thus we will focus on 
n > —1, where the danger of annihilations is in the past: 
T/H > 1 for 

T > Ta ~ io- 3 ( 3 +™)/( 1 +") GeV (m/GeV)" /(1+n) , (4) 

or 1 eV for n = 0. To ensure that early annihilations 
do not reduce CDM particles to negligible numbers, they 
must be protected against annihilation in the early Uni- 
verse. We suggest two mechanisms; doubtless, there are 
other possibilities. 

First, CDM particles could be produced late (T < Ta) 
by the decays of another massive particle. Note that this 
requires a long lifetime, r > t(Tj\) ~ 10 5 yrs, and the 
mass difference between the two particles should be small 
enough to ensure that the relativistic decay products do 
not make the Universe radiation dominated. 

The second way of avoiding the early-Universe anni- 
hilation catastrophe is to make the mass of the anni- 
hilation product be dynamical. For example, a phase 
transition that takes place at T < Ta could change anni- 
hilation from being kinematically impossible to possible 
if the mass of the annihilation product dropped below 
threshold after the phase transition (or if the mass of the 
CDM particle rose above threshold). A variation on this 
theme is coupling the annihilation produced particle to 
a scalar field, cf>, with ((f)) ^ 0. As (cf>) decreases, either 
quickly to zero as a result of a symmetry-restoring phase 
transition, or slowly as (4>) rolls to the minimum of its 
potential, the product particle's mass may drop below 
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threshold, opening up the new annihilation channel, at 
T <T A . 

Finally, the CDM annihilation products must not in- 
clude photons because their 7-ray flux would far exceed 
observational limits. For example for 1 GeV CDM par- 
ticles, the flux would be around 10 s cm -2 sr _1 s -1 , some 
ten orders of magnitude above the observed diffuse 7-ray 
flux at 1 GeV. 

Observational Constraints. We henceforth restrict our- 
selves to n = 0. The contribution of annihilations to the 
evolution of the density profile is given by 



d[p(r)/ PA ] /dt 



(5) 



where p A = m/(o~Ato) and to is the age of the Universe 
today. Assuming the initial || profile to be NFW, the 
resulting density profile is 



P(r) = Pi 



[x(l 



Ps/pc 



(6) 



where x = r/r s , and a core of constant density, p corc = 
Pa, is clearly evident. However, the mass loss due to 
annihilations results in adiabatic expansion of the core, 
such that the quantity M(r)r is left invariant JtJ. This 
expansion results in a lower core density and one can 
estimate that the ratio p CO re/pA ranges from about 0.1 
(dwarf galaxies) to about 0.3 (clusters). We have verified 
this by more detailed numerical work which allows us to 
determine p CO ic/ PA as a function of halo mass. 

We now turn to the observable constraints on anni- 
hilating CDM. A robust prediction of the s-wave anni- 
hilation scenario is that the cores are more evident in 
smaller mass halos, as can be seen in Fig. ||. So we first 
turn to the galactic satellites in the Milky Way group [Q , 
of which there are 11 known. 




FIG. 2. The solid curve is the NFW density profile. An- 
nihilation-modified profiles are labeled by the virial mass of 
the halo: 10 16 M© (cluster) and 10 12 M© (galaxy). 

For a 10 s M galactic satellite, the core radius pro- 
duced by annihilations is about 1 kpc, which is about 
the same as the cut-off radius induced by tidal forces. 
Most of the galactic satellites have large velocity disper- 
sions (~ 10 km/s) for their stellar content, which suggests 



that they are CDM dominated fTcj] . If so, their internal 
velocity dispersions indicate that p co t e — 0(1 GeV/ cm 3 ) 
& 

We also looked at dwarf spiral galaxies and LSBs. One 
must use these with caution since van den Bosch et at 
]l2[ have recently claimed that most of the Hi rotation 
curve data do not have sufficient spatial resolution to put 
meaningful constraints on the halo cusps. They do iden- 
tify three nearby galaxies which have sufficient spatial 
resolution - NGC 247, DDO 154 and NGC 3109 [H. 



van den Bosch et al. |12| find that 0.55 < a < 1.26 for 
the LSB (NGC 247), and a < 0.5 for the two dwarfs, at 
the 99.73% confidence level which at face value, argues 
for soft cores in low-mass systems. The annihilation sce- 
nario naturally explains this since the cores are more evi- 
dent in low-mass systems (see Fig. ||). However, it should 
be noted @ that the error bars on the rotation veloc- 
ity data are probably not a complete description of the 
total uncertainty and that a critical reevaluation might 
lead to a less stringent bound on a, thus alleviating the 
discrepancy between the observed dwarf rotation curves 
and CDM predictions. 

To estimate the cross-section required to achieve con- 
sistency with observations, we fit to the two dwarf galax- 
ies identified above with the halo profile in Eq. ^, a thin 
stellar disk and the observed gas. We have included the 
effect of finite resolution. We find that p A — 0.2 Mq/ pc 3 
results in a good fit to both (see Fig. ||). In both cases, 
the outer parts of the halo (determined by p s and r s ) are 
consistent with NFW theory. 
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FIG. 3. Rotation curve fits with p A ^ 0.2 M & / pc 3 . 

To see if structures on the largest scales are consistent 
with the annihilation scenario, we turn to strong gravita- 
tional lensing of background galaxies by clusters. Tyson 
et al. model the mass distribution in the cluster CL 
0024, which produces multiple distorted images of a back- 
ground galaxy, and find evidence for a compact soft core 
(of about 35/i -1 kpc) in the projected density. A similar 
value for the core radius was inferred earlier by Smail et 
al. Jl5| for CL 0024 and other clusters. X-ray studies 
(Bohringer et al. Jl5|] ) of CL 0024 are also in agreement 
with the above results. However, we urge caution in in- 
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terpreting these results since the evidence for soft cores 
in clusters is largely based on just one cluster (CL 0024). 

A value of pa = O.2M /pc 3 would produce a core 
density of about 0.06 AT Q /pc 3 in a cluster-sized object. 
We find that this core density is consistent with the sur- 
face density reconstruction of CL 0024 by Tyson et al. 
JL5| . The implied CDM mass within the arc radius (of 
107 /i -1 kpc) is in agreement with the quoted value of 
about 1.66 x 10 14 /i _1 M for the total mass within the 
arc radius (l(| . 

Discussion. The s-wave annihilation scenario with a 
cross-section of (a\v\) = 10~ 29 (m/ GeV) cm 2 produces 
a core density of about 1 GeV/ cm 3 over widely differ- 
ent scales. Intriguingly, this seems to be consistent with 
observations. 

Apart from cuspy cores |]l7f , simulations of non- 
interacting CDM also predict a much larger number of 
sub-halos for a galactic size halo than the observed num- 
ber of galactic satellites ]ll| . Certainly, the s-wave an- 
nihilation scenario has a dramatic effect on the small- 
est halos, and this could contribute to their destruction. 
However, further study is required to test this hypothesis. 

Another particle-physics solution in which CDM par- 
ticles have a large cross section for self interaction (a ~ 
10~ 24 cm 2 has been discussed jL9|. This possibility is be- 
ing tested by numerical simulations ]2p| ] . However, there 
are indications that self interactions lead to halos that 
are inconsistent with observations plj. The jury is still 
out. 

The requirements on a model for annihilating CDM are 
stringent, but by no means impossible ^2j. They point to 
a particle beyond those currently being considered, and 
therefore, to new physics. While it is possible that the 
solution to the CDM cusp problem will involve the inter- 
pretation of the observations or less exotic astrophysics, 
it is appealing to think that the properties of halo cores 
may teach us about the fundamental properties of the 
CDM particle. 
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